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The issue of averaging randomness is addressed, mostly in nuclear physics, but 
shortly also in QCD. The Feshbach approach, so successful in dealing with the 
continuum spectrum of the atomic nuclei { optical model), is extended to encompass 
bound states as well {shell model). Its relationship with the random-matrix theory 
is discussed and the bearing of the latter on QCD, especially in connection with 
the spectrum of the Dirac operator, is briefly touched upon. Finally the question 
of whether Feshbach's theory can cope with the averaging required by QCD is 
considered. 

1 The random-matrix theory and nuclear applications 

The random-matrix theory (RMT) is buih out of any stochastic modelHng of 
the matrix representation of a Hamihonian. It was introduced half-a-century 
ago by Wigner, who considered an ensemble of Hamihonians, different but 
sharing some symmetry, having independent random variables as elements 
and searched for those properties common to (nearly) all the members of the 
ensemble. For this purpose he introduced the probability of the occurrence 
of a given matrix H in the ensemble by multiplying a weighting function 
Pni3{H), — N being the Hilbert space dimension and (3 an index specifying 
the symmetry of the ensemble, — for the differentials of all the elements of 
the matrix representing H . 

Wigner focused on matrices with random elements having a gaussian 
distribution (Gaussian statistical hypothesis). The weighting function of a 
Hamiltonian with N quantum states is then 

P^,(iJ)^e-^^''«VA^^ (1) 

A being some A'^-independent constant. 

Concerning the symmetry, in the framework of the quantum Gaussian 
random-matrix theory (GRMT) one can define three different ensembles: 

i) the ensemble of the real symmetric matrices, describing time-reversal and 
rotationally invariant systems, referred to as GOE (Gaussian orthogonal 
ensemble) and corresponding to /? = 1; 
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ii) the ensemble of Hermitian matrices, describing systems violating time- 
reversal invariance, as, e. g., a nucleus in an external magnetic field, 
referred to as GUE (Gaussian unitary ensemble) and characterized by 



iii) the ensemble of the matrices that are linear combinations of the 2x2 unit 
matrix and the Pauli matrices, namely 



where Hnm and Hnm are real matrices, symmetric and antisymmetric, 
respectively. The matrices (|^) describe systems with half-integer spin, 
time-reversal, but not rotationally, invariant and the associated ensemble 
is referred to as GSE (Gaussian symplectic ensemble), wich corresponds 
to the choice /3 = 4. 

The elements of the matrices belonging to the ensembles GOE, GUE and GSE 
are real, complex and quaternion numbers, respectively. 

How to derive predictions on observables and how to compare these with 
the data in the XJRMT framework are questions recently reviewed by Weiden- 
mueller et ai. El These authors emphasize that predictions having universal 
validity, i. e. unrelated to the specific nature of the system under investigation, 
concern the fluctuations around mean values. Typical in this connection is the 
energy spectrum of a complex system with its mean level spacing A(e) and 
the related fluctuations: In the GRMT the latter are universal, the energies 
of the individual systems are not. 

More generally, are the local correlations among the eigenvalues (and the 
eigenvectors as well) that are accounted for by the GRMT, using A(e) as an 
input to fix the parameter N/ in (^) . As it is well-known, an amazing 
example in this connection is offered by the short-range fluctuations in the 
spectra of the atomic nuclei near the neutron emission threshold (^ 8 MeV). 
Indeed, given s — S/A(e) where S denotes the actual level spacing, the data 
of the nearest neighbour spacing distributions of the nuclear levels, p{s), are 
strikingly in accord with the GRMT predictions as shown in Fig. ^ 

It is of importance to realize that the strong repulsion displayed by the 
nuclear levels, when their excitation energy is sufficiently large and their dis- 
tance sufficiently small, in the GRMT is engrained in Eq. (|l]). Indeed, if the 
eigenvalues Ei, E2, ...En and the eigenvectors of the matrices of the ensem- 
ble are chosen as new independent variables, then (|l|) factorizes in two terms: 
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Figure 1. Nearest-neighbor spacing distribution versus the variable s defined in the text. 
The GRMT prediction (denoted GOE) is compared to the data (histogram). Also shown 
the result obtained with a Poisson distribution. (Taken from Ref. d). 



One only eigenvalue, the other only eigenvector dependent. The former reads 

N 

Pnp{Ei..En) ^ Y[\E^- EnfUdEi (3) 

m>n 1=1 

and clearly displays the levels repulsion in the Vandermonde determinant that 
multiplies dEi.AEN. 

Actually, there is still a long way to go from (|^) to p{s). It is remarkable 
that Wigner was able to cross it with his famous ansatz 

p{s) = ase-^''\ (4) 

which turned out to be very close to the exact predictions of the GRMT.i 
Note that the gaussian fall-off of (||) is not related to the gaussian nature of 
the ensemble, but directly arises from the Vandermonde determinant in (^. 
Applications of GRMT to QCD will be just touched upon in Sec. |. 

2 Feshbach's theory 

Feshbach's approach assumes the existence of randomness in complex systems 
whose internal structure is revealed by their excitation spectrum. On the other 
hand, in this view there are no systems completely ruled by randomness. The 
latter is defined through the relation 

{a{A)a{B))^{a{A)){a{B))^ (5) 
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Figure 2. (a) Elastic scattering of protons by ^■^Mo for energies close to the s-wave isabar 
analog resonance; (b) what one obtains by averaging the data of (a). (Taken from Ref. □). 



a being some physical observable of the system, dependent upon a parameter: 
The system is chaotic if the value of a in correspondence of the value A for 
the parameter is independent of the value a assumes for any other value of 
the parameter, jSay B. 

Feshbach's u view is that systems totally ruled by randomness do not exist 
or, equivalently, that coherent combinations of a can always be found such to 
allow the existence of correlations where now, in contrast with the discussion 
of Sec. ^, the correlations are dynamical, not statistical, they relate to order, 
not to disorder. 

The question then arises: How these dynamical correlations can be unrav- 
eled? The answer is: With appropriate averaging procedures. Before formally 
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defining the average, it is illuminating to experience how it works in a specific 
example: The elastic scattering of protons from the nucleus ^^Mo at vari- 
ous scattering angles and in the energy range of a few MeV. In Fig. |^a we 
show the data obtained with a poor energy resolution: They indeed appear 
to display the features of a random process (note that the values A and B of 
a generic parameter in here correspond to specific values of the energy). 
But in Fig. |b the same data (at 9 — 165°) are shown as obtained with a high 
energy resolution. What appears is astonishing: A perfect resonance typical 
of a highly correlated scattering process. 

The Heisenberg uncertainty principle AE At > fi/2, provides the key to 
understand. Poor energy resolution means large AE, hence small At or short 
time for the system to interact: The data appear random. On the other hand, 
small AE entails large At, hence a long time for the system to average out 
the fiuctuations: The results are then those of Fig. ||b. 

The necessity of averaging procedures to properly deal with randomness 
is thus clear. This program has been carried out by FeshbachO in his unified 
theory of nuclear reactions. We here extend his formalism to the bound states 
of many-body systems, focusing on the ground state, to complete Feshbach's 
program and to illustrate how it works. For this purpose we start by parti- 
tioning the Hilbert space associated to a many-body system in the P and Q 
sectors, respectively. 



2. 1 Theoretical framework 

As it is well-known, the splitting of the Hilbert space induced by the projection 
operators P and Q transforms the Schrocdingcr equation 

Hil) = Eij) (6) 

into the pair of coupled equations 

[E ^ Hpp){Pi;) ^ HpQiQi^) (7) 
{E ~ Hqq){Q4,) = Hqp{P^), (8) 

the meaning of the symbols being obvious. From the above the equation 
obeyed solely by (P-iA) is derived. It reads 

H(PV) = E(P^), (9) 

the P-space Hamiltonian being 

n^Hpp + HpQ I Hqp, (10) 

(4) 



paper: submitted to World Scientific on February 8, 2008 5 



with 

eQ = E-HQQ- Wqq (11) 

Wqq = Hqp \ HpQ. (12) 
il/ — J^pp 

It is of significance that although Eq. (|^) is not an eigenvalue equation, since 
the energy E also appears in the denominator of its right hand side, yet its 
solutions only occur for those values of E which are eigenvalues of as 
well; in other words, a one-to-one correspondence between the values of E 
allowed by (^) and (|^) exists (see later for a further discussion of this point). 
However, the solutions of (^ in correspondence to the various values of E are 
not orthogonal. 

Now, we assume the quantum deterministic aspect of nuclear dynamics 
to be embodied in the P-space, the chaotic one in the Q-space. Hence, the 
strategy of averaging over the latter follows, although, admittedly, some fuzzi- 
ness does affect this partitioning. To set up the averaging procedure we start 
by the recognition that the wave functions in the Q-space are rapidly varying 
functions of the energy E, viewed as a parameter classifying their ensemble. 

Accordingly, we average over this ensemble following the prescription 

/oo 
dEp{E,Eo,e)f{E), (13) 
-C30 

/ being a generic function to be averaged over the variable E with the dis- 
tribution p{E,EQ,e). The latter depends, beyond E, also upon the value Eq 
around which the average, — taken over a range of E essentially set by e, — 
is performed. A distribution convenient for our purposes is 

1 e'^'' 

p{E,Eo,e)^——-^ (14) 

ZTTl E — [Eq — e) — IT] 

which is indeed correctly normalized being 

dEp{E,E„,e)^l (15) 



(one should let 77 0+ after the integration has been performed). Note 
that Eq. (|lj) extends in some sense the Lorentz distribution of the optical 
modelQ to the situation of a zero width state. Hence the present formalism is 
especially suited to deal with ground states, which are of course stable: We 
shall accordingly focus mainly on the latter in the following- 
Now, in the Q-space the wave functions are found to baa'Q'S 

(g^) = —Hqp^o, (16) 
eg 
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"00 being an auxiliary function that in the end disappears from the formahsm. 
By averaging Eq. according to the prescriptions (^3|) and ( [1^ ) , one then 
finds that the averaged wave function of the nuclear ground state in the P- 
space (here denoted by the angle brackets) obeys the equation 

n{P^)^Eo{Pij). (17) 

In (|l^) Eq is the mean field energy and 

'H = Hpp + HpQ- — Hqp (18) 

is the mean field Hamiltonian. This can be further elaborated since the sin- 
gularities of the operator l/eg lie in the Imi<^ < half-planeO. Accordingly, 
one gets 

= 2^ '^^ E-{E^-e)-irjE- Hqq - Wqq{E) ^^^^ 
1 1 



Eq - e - Hqq - Wqq{E — Eq - e) Eq~ e - Hqq - Wqq{E = Eq) 

the last passage holding if the energy dependence of the operator Wqq is 
mild and if the parameter e is not too large (it should be not too small either, 
otherwise the energy averaging procedure becomes meaningless). 

The insertion of (^9|) into (^^ leads then to the following useful alternative 
expression for the mean field Hamiltonian 

n = Hpp + VpQVQP , (20) 

H/Q — e — ii 

where the energy dependent operators 



^PQ-gpQj /°7'/ (21) 
V ^0 - e - J^QQ 



and 



represent the residual effective NNJnteraction. The usefulness of the Eqs. (|20|), 
( pT| ) and (|2^ ) was realized in Ref. EJ, where it was noticed that with their help 
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the pair of equations (2.1) can be recast, as far as (PV') is concerned, into the 
form 

(i?-7Y)(PV) = ^pq(QV') (23) 
{E-Hqq){Q^) = Vqp{P^), (24) 

which is suitable for expressing the mean field fluctuations (the "error" ) . 

Indeed, by exploiting the spectral decomposition of the operator {E — 
7Y)~^ in terms of the eigenfunctions (/)„ of the mean field Hamiltonian 7Y,^ 
one gets from Eq. (p3) 



which, upon left multiplication by ((/>o|, yields 

(0o|P^) = M^^. (26) 

In the second term on the right hand side of Eq. ([25|), the prime stands for 
the omission of the n = term in the spectral decomposition. 
Next, the insertion of Eq. ( p5[ ) into ( p^ leads to 

m)^ ^ ^, VQp\c^o){um, (27) 

E - riQQ 

where the operator 

hQQ = Hqq + Vqp (;^-^) (28) 

has been introduced. Finally, by combining (|2^) and (p7|), one arrives at the 
equation 

E~E^^{MVpQ ^ \ ^Qp|0o), (29) 
E - riQQ 

which is the basis for computing the mean field energy error (or the fluctua- 
tions of the energy associated with randomness). 

Although Eq. (|2|) is valid for any choice of the projectors P and Q, its use 
is in fact appropriate when the P-space is one-dimensional, as it was indeed 



"We omit here the difficult proof concerning the normalization and orthogonalization of the 
eigenfunctions \4>n) of W. 



paper: submitted to World Scientific on February 8, 2008 8 



the case in Ref. where this choice was made for sake of simphcity. For a 
two-dimensional P-space, as we shall later see, one should rather single out 
two, rather then one, terms in the spectral decomposition of the operator 
(E-H)-^ on the right hand sid^oLEq. (|5|). 

We refer the reader to Refs. auB for a discussion on how the average of 
the square of Eq. ( |2^ ) is actually computed (the average of ( p9| ) of course 
should vanish) and on how the complexity expansion is organized. Here, we 
confine ourselves to the leading term of this fast converging expansion. 



2.2 Energy averaging 

To understand better the significance of the energy averaging distribution (Q) 
we show how it works in the simple cases of a bi-dimensional (A) and of a 
tri-dimensional (B) Hilbert space. 

A. Bi-dimensional Hilbert space 
Let Ixi) and \x2) be the two normalized states spanning the space. Here 
the only possible choice for the projectors clearly is 

P=|Xi)(xi| and Q=|X2>(X2|. (30) 

Then, by expanding the operator 1/(P — Hqq), Eq. (|^) can be recast as 
follows 



i?-|Xi)«ii(Xi|-|Xi)«i2(X2|;^£(^)"(|X2)(X2|)"|X2)aj2(Xi 



\P^) = 0, 
(31) 



which, upon multiplying from the left by (xi| and exploiting the idempotency 
of 1x2) (X2 1, simplifies to 

|2 1 



E - ail 



ai2 



(XilPV) =0, 



(32) 



E - a22_ 
where the shorthand notations 

On = (xi|-H"lxi), a22 ^ {X2\H\X2) and ai2 = (xi|i?|x2) (33) 
have been introduced. This equation is trivially solved yielding the eigenvalues 

ail + 022 ± V («ii " «22)^ + 4|ai2p , (34) 



which coincide with those of H. It helps notice that the eigenvalues (^) are 
also found as intersections of the hyperbola 

E^aii + ^^ (35) 



LO — 022 
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Figure 3. The eigenvalues of a bi-dimensional Hilbert space. The matrix elements of the 
Hamiltonian arc taken to be an = —2, 022 = 2 and ai2 = 3, in arbitrary units. The 
exact eigenvalues (squares) and the ones of the energy averaged Hamiltonian Ti (circles) 
are shown. They correspond to the intersections with the straight line E = lu + e: The 
former with e = 0, the latter with e = 7. The stability of the lowest eigenvalue and the 
upward shift of the highest one are clearly apparent. 



with the straight line E ~ uj. 

Also the energy averaged Hamiltonian (^0|) can be expressed in the basis 
spanned by xi 8'iid X2 and one gets the mean field equation 

|2 



E - an 



ai2 



i? - e - 022 

The latter is again trivially solved yielding 
1 



(xi|(PV)> =0. 



E±^ - 
^ 2 



an + 022 + e ± \J (an - 022 - e)^ + 4|ai2p 



(36) 
(37) 

51) (with E 



which now corresponds to the intersections of the hyperbola 
replacing E) with the new straight line E = oj + e. 

From Fig. ^, where the solutions E± and E± are graphically displayed, it 
clearly appears that, while 5^ the solution E+ is much larger than 
the more so the greater e is. It is thus clear that the averaging distribution 
(p^), while mildly affecting the eigenvalue of H lying in the P-space, drives 
away the one lying in the Q-space. 

B. Tri-dimensional Hilbert space 

The space is spanned by the normalized states |xi): IX2) and Ixs)- Now, 
two choices are possible for the projectors, namely 

IXi)(Xi| + IX2)(X2| (38) 
Q=|X3)(X3| (39) 
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CO 



Figure 4. The eigenvalues of a tri-dimensional Hilbert space: The case of a one-dimensional 
P-space. The following matrix elements of the Hamiltonian are taken: an = —4, 022 = —2, 
a33 = 3, ai2 = 2, 013 = —5 and a23 = 3.5, in arbitrary units. The exact eigenvalues 
(squares) and the ones of the energy averaged Hamiltonian fL (circles) are shown. They 
correspond to the intersections with the straight line E = u) + t: The former with e = 0, the 
latter with e = 10. The stability of the P-space eigenvalue and the upward shift of those 
belonging to the Q-space are clearly apparent. 
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Figure 5. The eigenvalues of a tri-dimensional Hilbert space: The case of a hi-dimensional 
P-space. The matrix elements of the Hamiltonian are taken as in Fig. 4 The exact 
eigenvalues (squares) correspond to the intersections of the straight line E = lu with the 
continuous curve. The eigenvalues of the energy averaged Hamiltonian Ti. correspond to the 
intersections of the straight line E = ui + e, with e = 10, with the dashed curves. Note the 
dependence upon e. The stability of the P-space eigenvalue and the upward shift of the one 
belonging to the Q-space are again clearly apparent. 
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and 



P^\Xi){xi\ (40) 

Q= |X2>(X2| + |X3)(X3|. (41) 



In both cases, Eq. (H) can be recast as follows 

{E-an){E-a22){E-a33)-\ai2\^{E~a33)-\ai3WE-a22)-\a23\^{E-an) = 0, 
which is the cubic equation yielding the exact eigenvalues. Note that Eq. 



(42) 



is easily obtained with the choice (2.2), because in this case the operator 



{E — Hqq) ^ is expanded in terms of the idempotent operator |x3)(x3|- Not 



so with the choice (2.2), because now {E — Hqq)~^ should be expanded in 
terms of the operator (4l|), which is not idempotent. Actually, the larger the 
powers of the latter are, the more cumbersome they become. Yet, also in this 
case it can be proved that Eq. ( ^2| ) holds valid. 

Let us now examine the solutions of Eq. (^): As in the previous bi- 
dimensional case it is convenient to display the solutions graphically. For 
the partition (pT^), one finds that they are given by the intersections of the 
e-independent curve 

- |ai2p(^-a33) , |ai3p(cj - 022) Qi2a1:3Q23 , 0^2013023 /.Q^ 
E=an-\ ^ / X \ :fr^-^ \ :frr-^ \ :fr^-^ — , (4^) 



where 



I?(w) = (tj - 022) (t^ - 033) - 1^23 (44) 



with the straight line E = cj + e, as displayed in Fig. ^, where the case e = 0, — 
which clearly provides the exact eigenvalues Ei of the Schroedinger equation, 
— is also shown. From the figure, it transparently appears that Eq = Eq, 
whereas Ei >> Ei and E2 >> E2, the latter inequalities being stronger 
when the parameter e is large. 



In the case of the partition (2.2), the solutions are given by the inter- 
sections of the curve obtained replacing 022 —* 022 — e in Eqs. ( p3| ) and (Q) 
with the straight line E = uj + e, as displayed in Fig. ^ We face here a new 
situation, since now not only the straight line, but also Eqs. ( ^ ) and (^^ are 
e-dependent. Yet, one again sees that for e = one recovers the eigenvalues 
Ei , whereas when e 7^ the intercepts occur for Eq = Eq and Ei = Ei, but for 
E2 >> E2- Hence, we conclude that the action of the averaging distribution 
(0) affects very little the eigenvalues belonging to the P-space, while pushing 
off the ones in the Q-space by an amount proportional to e. 
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2.3 The P -space 

Having defined the energy averaging procedure, to get further it is necessary 
to define the operators P and Q. For this purpose the natural candidates 
as building blocks of the P operator appear to be the eigenstates \(j)n) of the 
mean field Hamiltonian (|2^), defined by the equation 

■H\4>n) = E„\(bn). (45) 

Their finding requires, however, the solution of a difficult self-consistency 
problem. Hence, we make the simpler choice of viewing as building blocks 
of P the Hartree-Fock (HF) variational solutions, which are, e. g., trivial in 
nuclear matter, the system we shall consider. 

For sake of illustration we start with a one-dimensional P-space by setting 

^= Ixhf)(xhf|, (46) 

Ixhf) being the HF ground state wave function of nuclear matter (the Fermi 
sphere). Then, on the basis of (^), one derives the mean field equation 

^" = + Eo^l-E ' ^^^^ 

which relates the mean field {Eq), the HF (£^hf) and the true {E) energies 
per particle, and the equation for the statistical fiuctuation of the energy 

where 

E l(V'2p-2h|V|xHF)|', (49) 
2p-2h 

the bras (V'2p-2h| representing the two-particle-two-holes (2p-2h) states of 
nuclear matter, whose average energy per particle is £2- Thus all the quantities 
appearing in Eqs. ( |4^ ) and ( |4^ ) are per particle, including the parameter e 
and the residual effective interaction V (which are accordingly divided by the 
nuclear mass number A). 

Note also that Eq. (|4^) gives the "fluctuations" of the mean field energy in 
the first order of the complexity expansion, which is based on an organization 
of the Q-space in blocks of excited states of increasing complexity (see Fig. ^): 
Here, the contribution to the error only arises from the sector of the Q-space 
set up with the 2p-2h excitations. Moreover, although the states of the Q- 
space obey well-defined, coupled differential equations (see Ref. 0), we describe 
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4p-4h 



Figure 6. The partition of the Hilbert space of nuclear matter in sets of increasing com- 
plexity. The first box on the left defines the P-space, the second one embodies the simplest 
states in the Q-space and so on. 



them with the HF muhi-particle-muhi-hole solutions, an approximation not 
impairing the orthogonaUty constraint P ■ Q — 0. 

The sum in Eq. ( p9| ) is performed over the ensemble^ of the 2p-2h excited 
states lying in an appropriate energy range (in Refs. 13'EI taken to be fixed 
essentially by the^parameter e), whose number J\f2 can be computed using the 
Ericson's formulaM for the density of the spin J iV-particle-A^-hole nuclear 
states, namely 

where 

9 = 7;— and (j'^=T\-^, (51) 
I tp \ 0, Ji 

T being the nuclear moment of inertia, tp the Fermi energy, a = A/S MeV^"'^ 
and E the excitation energy of the system. 

Now, Eqs. ( |4^ ) and (p8|), owing to the double sign appearing in the latter, 
set up two systems, each one including two equations, in two unknowns. Two 
options are then possible in selecting the latter: One can choose either the 
ground state mean field and true energies per particle, i. e. Eq and — 
assuming the matrix elements of the residual effective interaction to be known, 
— or i?o and when E is experimentally known, — which is indeed the 
case in nuclear matter. 

We take the latter choice, requiring in addition the coincidence of the two 
i?o obtained by solving the two systems separately. Actually, and notably, 
both systems lead to the same formal expression for the mean field energy per 
particle, namely 

^0 = ^ {(i^HF + ^ + e) - V(^^HF - ^ - e)2 + 4/32} , (52) 
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which holds vaHd for e < E'hf — E, the right hand side of this inequahty being 
positive because of the variational principle. The above, when — > 0, yields 

Eo=E + e, (53) 

in accord with ( pT| ) and (p^), but also with (^8|), which, for V ^ 0, gives Eq = 
E: Indeed, when the residual effective interaction vanishes no fluctuations can 
occur and hence the parameter e should vanish as well. 

On the other hand, the two systems yield two different expressions for 
the sum of the matrix elements of the residual effective interaction squared, 
namely 

:Ei 



/3f = v(AV2^^ { [Ei{1 + /AV2) - (e + V^^2) - E^p 



(54) 



Ei{l + VW2) - (e + vW2e2) - E, 



HP 



MEi-Ehf) 



and 



Pi = VAy2^^^^ { [Eui-l + VW2) - (-£ + /^W2e2) + 



E, 



HF 



Eui-l + VW2) - (-£ + VW2e2) + £^1 



HF 



(55) 



AeiE^-E^p) 



both of which vanish in the limit e — > 0, in accord with the previous discussion. 

Formula (jSJ) , solution of the first system of equations (the one with the 
"+" sign on the right hand side of (^)), yields the value of the energy de- 
pendent quantity ( ^9| ) (here denoted by /3f) on the lower border of the energy 
band expressing the fluctuations of the ground state energy E (and thus en- 
compassing Eg). Formula (^5|), solution of the second system of equations 
(the one with the minus sign on the right hand side of (^)), provides instead 
( ^ ) (here denoted by f3l) on the upper border of the band (remember that 
E~e2< 0). 

Of course, of the energy per particle E we only know the experimental 
value, not the values on the borders of the band: As a consequence, we can 
only surmise the width W of the latter, thus providing two different inputs for 
the energy E appearing on the right hand side of (2.3), namely E^ — E + W/2 
and El = E — W/2. However, we can explore whether, for a given W, a value 
for the parameter e can be found (not too large, not too small) such to have 
the two mean field energies per particle to coincide. If this search succeeds, 
then an orientation on W (or, equivalently, on the size of the fluctuations of 
the ground state energy) can be gained. 
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Note that the framework above outhned holds because the same quantity 
0^ appears in both ( ^7| ) and (|4^). This |Occurrence stems from an approxima- 
tion whose vahdity is discussed in Ref. u. 

We now extend the formahsm by letting the projector P to encompass, 
beyond the ground state, the 2p-2h excitations of the HF variational scheme 
as well. Thus, instead of Eq. (El), we write 



P^\x^F){x^F\+y^\x^^){A (56) 



/3 



where the sum is meant to be extended to the whole set of 2p-2h HF excita- 
tions Ixhf)- 



With the choice ( p6[ ) the mean field Hamiltonian (20) is then defined 
and one can compute the mean field ground state energy per particle, Eq = 
{4>o\'H\4>o) , using for the ket |0o) the expression 

I'/'o)-.so|xhf)+E^2IXhf)- (57) 



7 



and accounting for the influence of the Q-space on the ground state mean 
energy per particle in first order of the complexity expansion, i. e. by setting 

Q = EIXhf)(Xh1pI: (58) 

7 

the sum running over the whole set of the HF 4p-4h excitations. In ( |57| ) sq 
and si are complex coefficients, fixed, in principle, by Eq. (^5|) and satisfying 
the normalization condition 

ko|' + El^2l' = l- (59) 

7 

After straightforward, but lengthy, algebra using (^^ and ( p8| ) one arrives at 
the following new mean field equation 

^0 = kopi^HF + E l^2p(XH'ipl^^lXH"F> + '^^l E 'I^XUfMx^^) + 
7 7 

' ^'^"'^Kx^fFl^lx^'^F)^ (60) 



Pi 

V being the bare NN potential. 

Notably, Eq. (|6^) turns out to formally coincide with (^). Indeed, the 
first three terms on the right hand side of Eq. (|6C| ) just yield the mean value of 
the original bare Hamiltonian H in the state (|57|). In the thermodynamic limit 
of nuclear matter out of these pieces only the HF energy survives. Indeed, the 
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correction to the HF energy per particle due to a finite number of particle-hole 
excitations vanishes in the thermodynamic limit. In other words, the 2p-2h 
admixture into (|5^) does not change the expectation value of the Hamiltonian. 
Hence, defining 

e = Y.\^'iWxt^^\v\x%)\\ (61) 

( |60| ) can be recast into the form 

^0 = ^HF + (62) 



whose similarity with Eq. (47) is transparent. 

The same will take place for any admixture of A^p-A^h HF excited states 
in Ixhf): Hence, in our framework different choices of the projection operator 
P lead to the same structure for the mean field equation for nuclear matter. 
This invariance does not hold in finite nuclei. 



The only, of course important, difference between ( [49[ ) and (61) relates 
to the residual interaction V , which in (^) induces transitions from 2p-2h to 
4p-4h states, rather than from the Fermi sphere to the 2p-2h states. 

Concerning the statistical fluctuation equation one can again use (|2^), 
with the state |0o) given now by Eq. (p7|). Then invoking the randomness of 
the phases oJJ.he. Q-space wave functions (RPA) and proceeding exactly as 
done in Refs. IZI'eI, one deduces the new fluctuation equation 

where £4 denotes the average energy per particle of the 4p-4h HF states. In 
(|63|), Mi represents the number of 4p-4h excitations contributing to the sum 
over the index [3 in (|6l|). 

Hence, the "formal invariance" , with respect to the choice for the projector 
P, holds for both the equations at the core of our statistical approach to 
nuclear matter in first order of the complexity expansion. 

Indeed, the inclusion of A^p-A^h states (with N > 2) into the P-space 
would merely imply the replacement, in (^3|), of J\fi with JVn+2 and, at the 
same time, to have defined in terms of the matrix elements of V between 
A^p-A^h and (A^-|-2)p-(A^ + 2)h states. In addition, one should of course insert 
in the energy denominator the average energy per particle of the (A^ + 2)p- 
(A^ + 2)h HF states. 

Therefore, the extension of the P-space rapidly leads to the vanishing of 
the fluctuations, owing to the very fast increase of the number A/jv- 
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However, as already mentioned, Eq. ( p9[ ) is not in general a good starting 
point to derive the fluctuation equation. Indeed, it selects out only one term 
in the spectral decomposition of the operator 1/{E — Ti) , which is appropriate 
for a one-dimensional P-space only. Hence, in place of (p5|), we rather write 

|P« . |*.)<*2™ +|^,<*2™ + (^)'V,,|W, ,64) 

with an obvious meaning of the double primed operator in the last term on 



the right hand side. Clearly, Eq. (64) does not follows directly from (pq), but 



it assumes that in the spectral decomposition of the {E — 'H)~^ operator only 
one prominent (collective) state \4>2) enters beyond the ground state \(t>o). 
Sticking to this model, instead of Eq. (p7|), we likewise write 



IW) = (2) VQp{\M{MP^^) + \(l^2){h\Pm, (65) 

being 



h^QQ = HQQ + VQp(^^) VpQ. (66) 



Then, by left multiplying ( |64| ) (with E = Eq) by {(I>q\ and using (|65|), we 
obtain 

(67) 



which generalizes Eq. ( |2£ 

In a perfectly analogous fashion, by left multiplying ( |6^ ) (with E ^ E2 
by (02 1 and using (p5h , we obtain 



(68) 

In the above, Eq and E2 stand for the first two exact eigenvalues of the 
Schroedinger equation; Eq and E2 for the corresponding quantities associated 
with Eq. (p5|). Eq. (^8|) shows that, in the present framework, all the energies 
of the P-space fluctuate. 

Now, the energy averaging of (|6^) and (|6^) vanishes by deflnition, but the 
energy averaging of their square, wJiich yields the "error", docs not. Hence, 



proceeding along the lines of Refs. un, we subtract on the right hand side of 
both equations their average values, square the expressions thus obtained and 
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make use of RPA, keeping of our expansion in the complexity of the Q-space 
states the first term only. Next, we exploit the structure of |02), which, like 
|(/)o)) must be normalized, orthogonal to |0o) smd of the form 



102) = >^cf 1X^1.) +co|xhf). (69) 
Finally, we arrive at the equation 



where, in addition to (pT]), the further definition 

e = T.^rc;\{xZ\v\x'^^)\' (71) 

/37 



has been introduced; moreover, we have set 

_ (02 I P^) 



(72) 



Likewise, for the energy of the 2p-2h state of the P-space one obtains the 
fluctuation equation 



where, naturally. 



£/2 — £4 V A/4 r 



v' = T.\'I\'\{xt^v\V\xll)\'. (74) 

/37 



We thus see that the statistical fluctuation equations (|70| ) and (|73| ) are actually 
coupled through the term (|7l|). 



Concerning the mean fleld equations, clearly with the projector (56) an 
equation should exist also for the energy of the 2p-2h state. It can be derived 
by computing E2 ~ (02 1^102) and, notably, it turns out to read 

^2 = 4V;g^ZT-^' 

E^p representing the HF energy per particle of the system in the 2p-2h excited 

(2) 

state. Since we split E'^p into a part associated with the HF ground state 
and a part associated with the 2p-2h excitation energies, both per particle, 
and since the latter vanishes in the thermodynamic limit, — as previously 
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noted in commenting Eq. (^0|), then E^p = E'hf , a relation we expect to be 
approximately fulfilled also in a heavy nucleus. 

We conclude from the above analysis that our approach leads to a set of 
mean field equations, one for each of the states lying in the P-space: These 
equations, unlike the fluctuation ones, are not coupled. 



2.4 Normalization and fluctuation of the P-space ground state wave 
function 

In the present framework the ground state spectroscopic factor S is the square 
root of the norm of \Pip), the system's ground state wave function projection in 
P-space. To find it we exploit the completeness of the normalized eigenstates 
of H. Hence we write 



M 



(76) 



Now, confining ourselves to set M = 1, then the Eq. ( |27| ) for IQtp) is warranted 
and we rewrite (u&i as follows 



= l-{^o\VpQ- 



1 



:VQp\^o)\{MPi')\' 



Moreover, when M = 1 then 

52 = |(0o|P^)p. 
Hence, by exploiting (|29|), Eq. ([ttI ) can be recast into the formB^ 



— {Eq - £'0) + _-^o ~ -^0 
dEo Eq — e — Eq 



Finally, employing (p2) the expression 



3 lE^F-EQ-e-2dp'^/dEQ 



2 yc^ 



HF 



Eq 



4/32 



Eq — e — Eq 



(77) 
(78) 
(79) 

(80) 



foUowsQ'u, where the energy derivative of the sum of the square moduls of 
the vacuum-2p-2h matrix elements of the effective interaction appears (its 
explicit expression is given in Ref. u). Note that ( |80| ) goes to one as (3'^ 0, 
as it should. 

If, however, the expression for the Q-space wave function aoDropriate for 
a two-dimensional P-space, namely (65), is used, then (see Ref. Ej for details) 
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one ends up with the expression 



d{E2 - Eq) 



"3^ 



E2 " Eq 

Eo — e — Eq 



dEn 



, (81) 

Eo-€-Eo ^ ^ 

which reduces to as it should, if {if)2\Pip) 0, i. e. for a one-dimensional 
P-space. In deducing the approximation 

(82) 



QQ 



Eo-h 



(2) 



has been made. 

Since, from Eq. (17: 



(83) 



then Eq. (|l]) can be recast as follows 
1 



5"= 1 



Eq — e — Eq 



1 



(E2 ~ Eq) 



,dE2 



1 + r2 \dEo dEo 



(84) 

which again reduces to (^ as r ^ 0. 

We now address the problem of the fluctuation of \Pt/j). For this scope, 
we focus on the ground state and, by combining Eqs. (p^), ( p^ ) and (24), we 
obtain 



{E - 7Y)[|PV) - l(PV^))] + (i? - Eo)\{Pi^)) = Fpgig^^) 



(85) 



(the angle brackets meaning energy averaging). 

Then, if use is made of the expression ( [27| ) for \Qtp) and of the spectral 
decomposition of the operator {E — H)'^ , one gets 



\p^) - |(p^)) = 



1 



S — En 
#o)(0o|VpQ 



~ ^QQ 

1 



E~h. 



QQ 



Vqp\^o){MP^) 



iE - E,) ^ (0n I (^^)>- 



(86) 



Now, since [(P?/")) cx |0o)i from the above finally it follows 



m-\{pi^)) = ii-\Mmr' 



E-n 



V, 



PQ 



E-k 



QQ 



-VQp\^o){<f>o\Pi^) 
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= i^olPi') (tt^) VpQ \ VqpIM, (87) 

which vanishes when P is given by (^^, since clearly \Pi^) does not fluctuate 
in a one-dimensional P-space. 

If, on the other hand, P is given by Eq. (p6[), then the above can be 
computed in first order of the complexity expansion, using ( ^7|) for |0o), ( |69|) 
for 1 02) and (|5^ ) for Q. One ends up with the expression (we set E = Eo to 
conform to previous notations) 

[ip^)-i(pv.))], = iMPi^)JM^Y.ix^^\y\xt^^)^^ixZ^^^^ 

Eq - E2 Eo - Ehf 

(88) 

which can be further simplified invoking the randomness of the phases of the 
wave functions in the Q-space and again introducing the average energy per 
particle £4 for the 4p-4h HF excited states. Then, with the help of Eq. (|7l|), 
the formula 

is derived. It gives the fiuctuations of the wave function in first order of the 
complexity expansion. _ _ 

Notice that in (89) the scalar product {(j)o\Pip), unlike in Refs. Q'EI, does 



not coincide with the spectroscopic factor S, as defined in ([76|), since Eq. ( p9D 
refers to a P-space with dimensions larger than one. Rather, it measures the 
amount of the true ground state wave function of the system embodied in the 
mean field state |(/io)- 



The relevance of Eq. (89) lies in the possibility it offers to assess the 
"error" affecting the ground state expectation value of operators associated 
to physical observables like, e. g., the magnetic moments of nuclei. 



3 A glance at QCD 

Lately the existence of randomness has been unraveled also at the scale of 
quarks, namely at the level of the physics ruled by QCD. Here the RMT has 
achieved outstanding successes in connection with the low-lying eigenvalues 
of the Dirac operator for massless fermions: 

^p = ^^ + g^^/l^ (90) 
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where g is the coupHng constant, the matrices A° the generators of the gauge 
group and the /l" the gauge fields. Indeed, the spectrum of (|90|), 

ip(l>n = A„(/)„, (91) 

has been computed for a SU(2) gauge theory (and for others as well) and from 
the A„'s the nearest-neighbor spacing distribution P{s) has been deduced. In 
turn, this has been compared to the prediction of the chRMT (chiral random 
matrix theory): An impressive agreement has been found J13 pointing to the 
universality of the fluctuations of the eigenvalues of the Dirac operator. 
The focus on the A„'s also stems from the anticommutator 

{*p,75} = 0, (92) 

whose vanishing in the chiral limit implies for the eigenvalues of (|^) the 
occurrence in symmetrical pairs (±A„) around A = 0. As a consequence the 
A„'s tend to accumulate around the origin and their average density 

p(A) = (^5(A-A„)), (93) 

n 

the average being taken on all the configurations of the gauge field, on the 
one hand displays the level repulsion already encountered at the level of the 
physics of the nucleus, — thus signalling universality, — and on the other 
plays a central role in the physics of the quark-gluon plasma. 

In fact, the density of the smallest eigenvalues of the Dirac operator_is 
directly related to the chiral condensate according to Banks-Casher formulaO: 

{qq) - -^p(0)/y, (94) 

V being the space-time volume. It is fascinating that such a fundamental 
aspect of nature as the breaking of the chiral symmetry appears to be ruled 
by chaos. 

In conclusion, we feci that RMT, in both its normal and chiral versions, 
is invaluable in disentangling the stochastic content of nuclear physics and 
QCD. In the former case, in the Feshbach's approach, stochasticity should be 
averaged out, thus leading to an effective Hamiltonian that has been proved 
most successful in interpreting the data. We argue that the same path might 
be worth following also in QCD. 
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